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FACILITY ABOARD SPACELAB AND A SUB-SATELLITE 

Peter L. Dyson* 

ABSTRACT 

The Spacelab program offers an opportunity to put a sophisticated Wave 
Injection Facility into orbit to study many aspects of space plasmas. Such a 
facility would ccver a very wide frequency range including MF and HF. In the 
Alouette and ISIS programs experiments were also done at these frequencies. 
Significant technical developments have been made since then, however, and it 
is now possible to carry out more extensive measurements on the properties of 
waves propagating in the space plasma. For instance, phase and Doppler shift 
measurements can he made to give more accurate measurements of echo time 
delay and the angle of arrival. This report examines various modes of opera- 
tion of a Wave Injection Facility in order to demonstrate the type of measure- 
ments which can be made. Since the Spacelab flights will Involve sub-satellites, 
some consideration is given to propagation between two vehicles both at rtF 
and VHF. 
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OPERATIONAL MODES FOR A WAVE INJECTION 
FACILITY ABOARD KPACKLAB AND A SUB-SATELLITE 

1. Introduction 

This report considers the operation of a Wave Injection Facility primarily 
from the |joint of view of meusuring the properties of electromagnetic signals in 
the MF and HF hands. The modes of operation suggested need not, however, be 
restricted to theBe types of sinful ami can he used in other frequency ranges 
and for other types of signals (such as electrostatic waves associated with plasma 
resonances). It should also be noted that the operating modes suggested here are 
not necessarily the only desirable modes at MF and HF. They are meant to be 
basic general purpose modes of operation which will yield a considerable amount 
of useful information. In fact, other modes of operation (or variations of the 
modes discussed here) will also be used for certain studies. Mow'ever, it is felt 
that the oj»erating modes presented are sufficient to illustrate the MF and HF 
capabilities of the proposed facility. 

It should be understood that the type of computer controlled facility dis- 
cussed here offers great flexibility in designing modes of operation and that new 
operating modes can be implemented at any time simply by writing a new control 
program, i.e., only changes in software are required. This type of computer con- 
trolled operation has proven successful with ground based ionosondes (e.g. Bibl 
and Reinisch, 1978; Devlin, 1978). The development of the Wave Injection Facility 
will be evolutionary i . nature and a flexible computer controlled facility is con- 
sistent with that philosophy. 
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A specific case of combining II F und VHF measurements to measure proper- 
ties of gravity waves is discussed. Again, this is an illustrative example and is 
not meant to exclude other modes of combined operation. 

Finally, in developing modes of operation one is usually concerned with the 
sequence in which certain measurements are made ra* h, 'r than the actual hard- 
ware equipment required to make the measurements. Con equently, no details 
will be given here of the equipment to be used. Ilather the emphasis will be on 
the measurement sequence and the methods of combining measurements to de- 
termine useful quantities. 

For convenience most of the discussion will lie in terms of a monostatic 
system but all of the concepts are appropriate to the bistutic situation. The 
electromagnetic wave studies at MF and HF will be primarily concerned with 
tran emitting pulses and determining the time delay, amplitude, phase and Dop- 
pler shift of any echoes received. As pointed out by Calvert (li)7(i), coherent 
detection of echoes preserves these signal properties and if the In phase and 
quadrature components of the coherently detected signal are digitized, further 
signal processing can be done by a computer or dedicated microprocessors. 
Details of actual receivers, coherent detectors and AGC circuits will not be 
given here since such devices arc routinely built for communications and radar 
applications. 

2. Amplitude Measurements 

After coherent detection, the in-phase and quadrature components of the 
signal are available for processing. Because of the large range of signal ampli- 
tudes which can occur it is unlikely that a detector with sufficient dynamic range 
could be easily designed. The standard approach is to adjust the receiver gain 
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Cfing an AGC system so that the input signal to the detector will be at an appro- 
priate level. Thus the measurement of signal amplitude requires the measure- 
ment of two quantities, the AGC level (receiver (tain) and the output level of the 
detector. It is unlikely that the AGC will be adjusted more than once for each 
transmitted pulse and it may be adjusted at a slower rate. To obtain an upper 
limit on the data rale from the experiment it will lx* assumed that the AGC level 
is measured at the rate at which K F pulses are transmitted. For digital proces- 
sing it is convenient to step the AGC by discrete amounts. 

Suppose that for each AGC setting, the output ol the AGC is divided linearly 
into 25<> levels so that it can he represented by an 8 bit word. The dynamic range 
at each gain setting is then 48 db (20 log, 0 256). If we assume an overall dynamic 
range of 200 db is required, the AGC control must cover 150 db. If the AGC level 
is changed when the detector output changes by a factor of 10 1 8 , the AGC will 
have 00 levels which requires 6 bits. The AGC is changed when the maximum 
signal arriplitudc changes by the relatively small amount of 10 1 8 (2.5 db) to en- 
sure that most of the full dynamic range of the detector (48 db) is used for each 
transmitted pulse. By choosing a factor which is a power of 10, the 8 bit words 
which represent the signal amplitude (the AGC level and the detector output) are 
easily combined to give the logarithm of the signal amplitude. 

3. Phase Measurements 

The phase may be obtained from the in-phase and quadrature components 
of the detected signal. If we present the complex detected signal by z » x ♦ jy 
where x and y are the in-phase and quadrature components, the phase of the 
signal is given by 

tan 0 = 
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The error In ; , 64 Is related to the errors, ft x and 6 y. In x and y, by 


i; s * i n 4 c os 4 



Now for the same error in x and y. (i.e. ix»i y) we have 


6 0 = s in t coi ; 



* sin <t> ( * ♦ rot / ) — 


(sin 


. . &x 

• CO* • ) — 

A 


where A is the amplitude of the signal. For a given *x/A, the maximum error 
in / occurs when 4 • 45° 

&4 » — 

■■■ A 

11 the phase is to be measured to within 10 , then *x/A £ .17. Now the 
output of the detector is quantized into 250 levels so the quantization error is 
1/512 of full scale. This quantization error is therefore less than 17% for all 
signals w ith amplitudes greater than l.J?” of full scale of the detector output. 
Consequently, for the criteria previously outlined for amplitude measurement, 
the phase of signals can be measured to better than 10° over a dynamic range 
of :is db (-20 log,,, .012 ) for each AGC setting. Thus the procedure outlined 
for amplitude measurement also provides acceptable phase measurements. 

II a 12 bit word were chosen for each amplitude, then the dynamic range at 
each gain setting would lie 72 db and the phase wou.d be measured to within 10° 
over a dynamic range of 45 db. 
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Ground bused measurements with un ionosonde have shown that the phase can 
be consistently measured to within 10° (Devlin et al., 1977). It ma> be worth not' 
in« that Devlin et al. measured the phase more directly, by using limiting detec- 
tion of echo signals. Such an approach could also be used for the Wave Injection Ex- 
periment; the amplitude and phase, however, still require two computer words for 
their representation. Final choice between the two possible methods of measure- 
ment depends on hardware design considerations. Doth may l>c considered "state 
oi the art" techniques and the final choice does not affect the considerations here. 

4. Group Path (Apparent Range) Measurements 

If pulses of length r sec. are transmitted, the corresponding receiver band- 
width is 1/r. In order to preserve the signal information after detection, each 
component of the detected signal must be digitized at a rate of 2/r. If the detector 
out|>ut is continuously digitized at this rate the group path will be measured to 
within ct/2. In the monostatic situation, the group path will be twice the apparent 
range to the echoing region. Hence the resolution in apparent range w ill be 0.25 
cr. Thus for a bandwidth of 10 kHz, pulses of 15 km in apparent range would be 
transmitted and the apparent range quantization would be 7.5 km so that the ap- 
parent range could be measured to ±3.75 km. 

Range Vernier . The above range resolution may be improved by combining 
measurements made on two adjacent frequencies. The basic definition of group path 
path. P' is (e.g. Budden, 1961) 

p '=17< fp > 
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where I* id the phuse path ut frequency f. Now the phase of an echo signal, ; is 
related to P by 


$ * 


2 11 * 

c 


If the phase is measured on two frequencies differing by f, then 


A/ 2rrA(fP) 4 

s — . I - - - — 

Af v Af C 


R' 


where IV is the apparent echo range. Now the measured phase difference etui 
only be between 0 and 2 but the actual phase difference may be greater than 2 - . 

Thus .A; « 2m ♦ ; where n is an integer and - / is the measured phase 

difference, and 


R* 


1 c <• 

7 A f + 4 I 


(4.1) 


Consequently, only the second term in (4.1) can be determined by the measured 
phase difference between two frequencies. The maximum value of the second 
term is 1/2 c/A f. If A f ■ 20 kHz, the maximum value is 7.5 km. Since the range 
can Is.* measured to the nearest 7.0 km by the method outlined previously, the 
first term in (1) can Is* determined. If the phase difference can be measured to 
within*- 20 , equation (4,1) can be used to determine IV to better than 0,5 
km. 

A comment on approximating the derivative ._*// • f by Vt/Af is appropriate. 
Assume that the phases / , and t . ♦ A; at the frequencies f, and fj ♦ Af, are 
known. We choose to calculate IV at the frequency, f mid-way between f ( and 
f, ♦ Af. 

i.e. f » f, ♦ A f/2 ■ fj ♦ if. 


G 


Then, using Taylor's theorem we may write 


* " f * ‘ n • *• - ' ♦ i (0) #i < ‘ * * - t ( 3 \< ‘ 1 >’ ♦ • • • • 

*,*t*-«t*th-** (^m 4(||f) <‘o»* | ($) <»*>*♦ ... 

Now, the error in eouating the L.H.S. to the first term on the II. II. 8. 1 b 

(ai, 1 = i i'*.' (ar) 1 
24 \df'L hc an 

t ( 

and this term will generally be quite small. 

The discussion of the Range Vernier technique has ignored the phase change 
between transmitted pulses which occurs due to the satellite motion, (i.e. the 
Doppler effect) an. I this poini is discussed laU r in the section on Combined Dop- 
pler and Group Delay derived from Phase Measurements. 


Range Window . The time between transmitted pulses defines a range window 
for the apparent range. Echoes arriving from greater apparent ranges will either 
appear at a shorter range as determined by the time difference between the most 
recently transmitted pulse and the echo, or will not be detected if the receiver 
frequency has been changed by more than the bandwidth of the receiver. For the 
proposed Spaeelab orbit, a range window of 500 km is appropriate in most in- 
stances. This corresponds to a transmitter pulse repetition rate (prf) of 300 
pulses per second (pps). In the equatorial region, long range echoes (up to 5000 
km or more) occur as a result of field-aligned propagation along spread F bub- 
bles and other irregularitie3. Such ranges can be detected uinambiguously by 
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using a p.r.f. of 30 p.p.s. It will be .shown later that such a p.r.f. is too low to 
muke tin* required Doppler measurements. Range ambiguities produced by 
higher p.r.f.'s can in* removed in two ways. Firstly, the Doppler shill may Ik* 
used since the longer delay echoes will generally have a different Doppler shift 
than the shorter delay vertical incidence echoes. However, this method requires 
some knowledge or assumption of the propagation paths involved and may not 
always reliably remove range window ambiguities. Secondly, coded |>ulses could 
he transmitted and the code varied from pulse to pulse and over a sufficient num- 
ber oi pulses so that no ambiguities remain. Each code would require its own 
detection system. This approach would not add significantly to the hardware of 
tlie system, but it would increase the amount of data to be digitized. Each addi- 
tional pulse code would require an extra in-phase and quadrature amplitude to 
be digitized. If a p.r.f. of 300 p.p.s. is used and echoes with ranges of, say, up 
to 3000 km are to be unambiguously measured, six pulse codes would be required. 

3. Doppler Measurements 

One standard method ol Doppler measurement is to pass the detected signal 
through a series of narrow band filters, the output of which correspond to the 
Doppler spectrum of the signal. While this method gives the instantaneous Dop- 
pler directly, there are two distinct disadvantages. The first is that the amount 
of data to be digitized is proportion; 1 to the number of filters and hence will in- 
crease rapidly if good Doppler resolution is required. The second is that the 
filter characteristics need to change with the transmitted frequency and this 
would complicate any hardware approach. (The reasons for requiring the filter 
characteristics to change with frequency will be apparent later). 
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An alternative approach is to measure the phase of an echo signal as a func- 
tion of time and to calculate the Doppler spectrum using the Fast Fourier Trans- 
form (FFT). This method has two main advantages. Firstly, it uses the phase 
information already available from the echo amplitude measurements and so does 
not increase either the hardware required or the number of measurements to be 
made. Secondly, at a given transmitter frequency, the frequency range and reso- 
lution of the Doppler spectra will be determined by the transmitter prf and the 
total time for which measurements are made. These parameters may be readily 
changed by software in the computer control system and so the required changes 
in Doppler resolution as the transmitted frequency changes can be more easily 
made. 

It may seem a disadvantage that this method combines measurements over 
some period of time and does not give an instantaneous Doppler as is obtained 
with hardware filter circuits. However, such filter circuits involve time con- 
stants, etc., and some signal integration is required to obtain appropriate fre- 
quency resolution and consequently the "instantaneous" values obtained are not 
based exclusively on the signal components received at a given time. 

Interpret at i on of Doppler 

The general formula for the Doppler shift in a refractive medium, is given 
by (Bennett, 1968) 

Af = --“ |J &co.ad. ♦!*.*]■} (5.1) 

where n is the refractive index 

p is a vector in the direction of the wave normal with magnitude, ^ , 
v is the vector velocity of a point on the ray 
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A ami B are the end points of the ray. 

Consider Initially an ionosphere which is plane stratified (but not necessarily 
in the horizontal plane) and which is not varying with time. For an orbiting ex- 
periment, the Doppler effect w ill be due to the relative motion of the orbiting 
vehicle and the ionosphere. For rays reflected back to the vehicle in the usual 
way (Figure 1), the Doppler shift is given by (e.g. Dyson, 1975) 



(5.2) 


2 f 

= — M . V COS 0 

C A 


(5.3) 


where im a is the refractive index at the vehicle, v is the vehicle velocity and 
the angle between the wave normal direction and the vehicle velocity. 

If the reflection point is not at the ^ * 0 level, but is due to an Irregularity 
at some other level, several return paths are possible (Figure 2). 

For paths 1 and 2 in Figure 2, the Doppler shift is still due entirely to the 
motion of the end point A and (5.2) still applies and can be used to obtain & x and 
2 . For rays outgoing along 1 and returning along 2, or vice versa, 

A f =7 ( p A + P A ) ’ y (5.4) 

t iA, -Aj 


= ~ (m Aj cos 6 X + m a , cos 0 a ) (5.5) 

and so equals the average Doppler shift for the signals returned by ray paths 1 
and 2 . If the ionosphere has a velocity v^, then v must be replaced by v - in 
*5.2) and (5.3). However, the corrections will generally be minor since |vj >> I v , • 
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For the more general ease in which the background ionosphere is not hori- 
zontally stratified, equations (5.2) and (5.3) still apply if d^/dt is approximately 
constant along the ray path (Dyson, 1975) and this situation is likely to hold for 
large scale features, such as TID's, which produce refractive effects. 

Hence, monostatic measurements of Doppler shift from an orbiting vehicle 
may be regarded as measurements of angle of arrival provided the refractive 
index at the satellite and the vehicle velocity are known. The earth's magnetic 
field causes the refractive index to depend on the wave normal direction and 
also causes the ray direction to differ from the wave normal direction. These 
relationships are well known and may be accounted for. In many cases, however, 
these complications may be neglected provided the frequencies under considera- 
tion are not near the frequencies reflected at the orbiting vehicle. The plasma 
frequency and gytofrequenev at the vehicle are measured bv the WIF and may 
be considered sufficient information to determine m a . 

The ability to measure angle of arrival is very important. It will enable 
the nature of ionospheric structures to be determined more precisely. An ob- 
vious example is ducted field aligned propagation which occurs generally near 
the magnetic equator (Muldrew, 1963) and can occur in association with equa- 
torial plasma bubbles (Dyson and Benson, 1979). Angle of arrival information 
would give confirmation that the propagation is in fact field-aligned and would 
allow the trapping mechanism to be studied in more detail. The question of 
whether extra traces usually attributed to conjugate ducting are always due to 
ducting or are sometimes due to simple refraction (I.obb and Titheridge, 1977) 
would also be resolved. 
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Kquation (5.2) may be re-written as 


Af 

f 



v cos 8 . 


Assuming a vehicle velocity of 8 km/sec, and >. A ~ 1 for frequencies greater 
than 1.5 times the plasma frequency at the vehicle, a Doppler resolution ■ f/f 
of 0.0 x io* ft gives an angular resolution about ' ■ 0° of ~10°. Near^ ■ 00° 
(vertical incidence) the angular resolution is -*1°. This resolution in Doppler 
corresponds to a change in Doppler of 1.8 Hz at 3 MHz and 0 Hz at 10 MHz. To 
achieve these resolutions measurements would have to be made lor at least 0.50 
seconds at 3 MHz and 0.17 seconds at 10 MHz. 

The FFT method will be used to calculate the Doppler from the phase meas- 
urements vs. time. The highest Doppler shifts will be of the order of 150 Hz at 
3 MHz and 550 Hz at 10 MHz. To prevent aliasing, the phase must be measured 
at a rate of 300 samples per second at 3 MHz and 1100 samples per ~eeond at 
10 MHz. As discussed earlier, a p.r.f. of 300 provides an appropriate range 
window (500 km), but the smaller range window associated with significantly 
higher p.r.f.'s will be generally unacceptable. It is possible to use lower p.r.f. 's 
and remove the Doppler aliasing by transmitting first at one p.r.f. and then at 
another. For example, suppose a measured sequence is made at a particular 
p.r.f. There will be an ambiguity for a Doppler frequency greater than half the 
p.r.f. If a second measurement sequence is obtained at a p.r.f. of 0.67 of the 
first p.r.f., the ambiguity can be removed since the ambiguous frequencies of 
the first spectra will appear as distinguishable apparent frequencies in the 
second spectra. 
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Suppose that the aliasing irequcncy for the first p.r.f. is f B . An apparent 
frequency f E) may be any frequency f which satisfies the relation f ■ Znf, * f E , , 
n ■ 0, 1, 2, . . . For the second p.r.f., the relationship may be written f ■ 
2n(.67)f a i f„ 2 , n»0, 1, 2, . . . It is apparent thut for a given f, f, ( ■ f, 2 only when 
n ■ 0. If a Doppler frequency f produces spectral peaks at f M( and f #J in the two 
spectra, then 

2nf. if, lS 2n(0.67)f. if„ 

and 

_ 

0.66 f, 

For example, suppose a peak occurs at 0.5C f„ in the first spectra and at 0.16 f, 
in the second, 


n = 


-0 34 
0.66 


or 1 


Since n must be a positive integer, n * 1 and f ■ 1.5 f„. If the higher p.r.f. is 
chosen to be 300 p.p.s., then for frequencies up to 10 MHz, nl 4. Of course 
if many peaks occur in the spectra some ambiguity in interpretation may remain. 


Simple Doppler Measurement 

If two consecutive pulses are transmitted on the same frequency and the 
phase difference measured, the Doppler shift may be calculated directly. This 
will give a reliable Doppler measurement for discrete echoes with a large sig- 
nal to noise ratio. When noise is present the method may still be used by trans- 
mitting several pulses and determining the average Doppler shift. When more 
than one echo is present at the same apparent range this method will give some 
weighted average depending on the amplitudes and Doppler shifts of the individual 


13 



echoes. In this case the lull Doppler spectrum must be measured In order to 
resolve the echoes. However, If one echo has significantly greater amplitude 
than other echoes the simple Doppler measurement will be a reasonable approx- 
imation. The Bimple Doppler measurement will be a very useful technique be- 
cause well defined discrete echo traces often occur and the technique is ideally 
suited for use with a scept frequency mode of operation. 

The Doppler shift and phase are related through the follow ing expressions 
<P * 1): 


2f f dP 

. f = — v . cos b - — —— 
c c (Jt 


4 - 2tt| P c 


Now 


■V _ 2rrf AP 
At ~ c At 


4nf 

c 


cos 6 


Now, for a p.r.f. of 300 p.p.s., and a vehicle velocity of 8 km/sec, the phase 
change between pulses is 

A<*>~ . 36 * 10* 6 f cos 6. 

The measured phase difference will be in the range of ±v , so that 

A<£ = 2i w + b<$; 

where / is the measured phase between *r;. Consequently, the angle of arrival 
may be aliased. The maximum apparent angle, c m will be given by 

I cos d | = — = 1 for f - 2.8 MHz. 

m 3 . 6 f 

Thus there will be no aliasing in the angle of arrival for frequencies less than 
3 MHz. (This result is expected in view of the discussion in the previous sec- 
tion). At 10 MHz, | cost m | = 1/3.6 and & n « 74° so that 90° * 6 . pp , r ,nt ~ 74 °- 
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For discrete echo traces, the apparent ungle of arrival will vary as a func- 
tion of the transmitter frequency In a predictable way enabling ambiguities to be 
removed. 

If the phase difference can be measured to within 20°, the resolution in angle 
of arrival, &<?, will be given by 

cos b0 | a 3 . 4 x 10V.36 f 
^•0.1 nt 10 MHz 

This corresponds to~6° at'.' * 90° (vertical incidence) and ~16° about 9 ■ o. At 
higher frequencies the resolution increases but aliasing becomes more of a prob- 
lem, Figure 3 shows the variation of A; with frequency and angle of arrival for 
a p.r.f. of 300 p.p.s. 

In Figure 4 a situation is shown which can occur near the equator. The 
Spacelab is in the topside ionosphere aad is within a field-aligned duct which 
produces conjugate ducting. Vertical incidence echoes are also present, as well 
as backscattering echoes from a patch of field-aligned irregularities below the 
vehicle. The echo patterns which would be detected using a large range window 
are shown. In Figure 5 the type of record obtained with a 300 km range window 
is shown. The echo intensity has been coded with the simple Doppler measure- 
ments (i.e. A<£) [in the manner first used by Findlay (1951)], to produce a fringe 
pattern in each of the discrete echoes. (Because Atf changes slowly with fre- 
quency [as seen in Figure 3], it has been multiplied by about 5 to produce the 
fringe pattern shown). The vertical incidence echo and the L+S ducted echo 
(see Figure 4), have zero Doppler shift and the fringes are parallel to the echo 
contour. The S and 2S+I, ducted echoes have a positive Doppler and the fringes 
move upward as a function of frequency with respect to the echo contour. For 
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thi* L ducted echo, the Doppler in negative ;ind ho the fringes slope downward. 

The backscattering from the spread F patch will generally consist of overlapping 
echoes and coherent echo traces will be difficult to distinguish. However, the 
closest range edge of the spread F echo patch will correspond to echoes from 
the edge of the spread F irregularity region and may show fringes of negative 
slo|H.' corresponding to propagation |>er|>endiculur to the magnetic r .eld direction. 

Combined Doppler and Group Delay Derived from Phase Measurements 

In a previous section it was pointed out that a range vernier technique could 
be used by measuring the phase difference on two frequencies. Since these phase 
measurements would be made sequential^ , the phase difference would include a 
Doppler shift contribution. This contribution may be removed by transmitting a 
sequence of three pulses, the first two on one frequency to measure the Doppler 
and the third at the second frequency, for the range vernier technique. 

We may write: 

dl _ dd '4<f. ^i_f 
dt <3t df dt 

where d^/dt is the measured rate of change of phase between two consecutive 
indses, •//< it is the Doppler contribution, df/dt is the sweep rate of the trans- 
mitter, and cV/ f is the contribution associated with the path of propagation. 

If / n , ; 12 , and 4 2) are the measured phases of the three consecutive pulses 
transmitted at intervals of A t and at frequencies f |t f, and f, ♦ Af, respectively, 
then 

(d£\ m -*n ♦ 2nv 

\ dt j , ~ at At 
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and 


/d£\ _ t *4 df . *ai "*ia ♦ 2mv 
\dt / 2 't Tf dt “ At 

Hence 

-(*1) 

of dt \dty, \dt /, 

* 21 • 2< *» * * II * 2 < ir " n )" 

__ 

<t> u - 2<£ 12 ♦ <*,, ♦ 2(m - n) w 
df ~ 

Hence both dtf/ct and ^f./~of may be determined with the usual 2t\-n ambiguity, 
the effect of which has already been discussed. 

Since three phase measurements are involved the errors will be greater 
than before, but with the values adopted in the section on "Range Vernier" the 
apparent range would be determined to within 1 km. 

6. Pulse Synthesis 

Although the Range Vernier technique allows the group pa*h to be readily 
measured to within 1 km, the ability to separate closely spaced echoes is de- 
termined by the length of the transmitted pulse, assumed here to be 15 km. If 
better resolution of echoes is required, pulse compression techniques must be 
used if the receiver bandwidth is not increased. One simple method has been 
used successfully with ground based ionosondes (Hammer and Bourne, 1976; 

Devlin et al, , 1977). The ionosonde is operated in the usual way in that a sequence 
of pulses is transmitted consecutively over an appropriately limited frequency 
range. The amplitude and phase of any echoes are measured. The echoes associ- 
ated with a particular echo trace (e.g. one hop vertical incidence F-region echoes) 
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may then be regarded a* *he "frequency components" of an echo signal due to a 
synthesized pulse whosi *quency spectrum is tin sequence of transmitted |>ulses. 
The delayed synthesized echo is then obtained by Fourier transformation of the 
amplitude and phase of the "frequency components." Quite narrow pulses can be 
synthesized and echoes separated by 1 km In range readily resolved (Devlin et al., 
l‘J77). The method is discussed in detail by Devlin et al. (1U77). It is an indirect 
form of pulse synthesis but is discussed hero because, for the type of computer 
controlled sy stem envisaged, its implementation depends only on appropriate 
software. It could also conceivably be used in later analy sis of certain swept 
frequency modes of o|>erution not specifically designed with this technique in 
mind. If pulse synthesis is required as a primary operation mode, however, 
more conventional methods should probably be used since the technique described 
here is relatively time consuming both from the |>oint of view of measurement 
time and the amount of analysis required to obtain good echo resolution. 

7. Location of Echoing Regions 

Vertical incidence echo traces may be inverted to obtain the electron density 
distribution with height. The procedures may also be applied at other angles of 
incidence and this has been done by Dyson and Benson <1‘J78), for example, to 
determine the electron density distribution along magnetic field-aligned bubbles. 

In the absence of angle cf arrival measurements, this type of analysis is usually 
not attempted when it is obvious that the angle of arrival is changing significantly 
with frequency, e.g. when TID's cause echo traces which exhibit "kinks" or 
splitting. Angle of arrival information w-ill be available in the Spacelab experi- 
ment and such techniques, combined with ray tracing, wdll enable echo traces 
to be inverted to obtain the irregularity structure. 
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An approximate method o( locution Irregularity regions distant from Space- 
lab will be discussed in the section on Fixed Frequency operation 

8. Operation Modes 

Swept Frequency Mode . This mode will be a basic exploratory mode aimed 
at getting a reasonably complete picture of the ionosphere in a short time. Thus 
both range and angle of arrival (Doppler) information are required, the lutter 
be inn obtained by the Simple Doppler Method. The principles of o|>eratlon will 
lie those outlined in section 5, i.e., a group of three pulses, two at the same fre- 
quency f and u third at a slightly higher frequency, f ♦ Af. A frequency inter\al 
Af of 20 kHz will give ranges within 1 km using the Range Vernier technique. 

Suppose 3 sets of pulses (i.e. a total of nine pulses) are transmitted to im- 
prove operation in the presence ol noise, and 'hat f is then increased by 25 kHz 
if below 1 MHz and by 100 kHz If above 1 MHz. If a p.r.f. of 300 p.p.s. is used, 
the frequency range from 500 kHz to 30 MHz will be swept in 4.5 Beconds. Gen- 
erally, ionospheric echoes will occur ovtr a much smaller frequency range than 
this and it will not lie necessary to cover the whole range of 500 kHz to 30 MHz. 
Consequently, very rapid soundings cun be made and the resolution of the meas- 
urements along the orbital path would be better than 35 km. Above 1 MHz, range 
measurements at every 100 kHz in frequency are sufficient for calculating the 
electron density distribution by inversion techniques. 

"Fixed Fre quency" Measurements . In this mode measurements are made 
at a small number of frequencies, say up to 10. This mode can be considered 
to be a variation of the previous mode since a certain frequency range must be 
cycled threugh using discrete steps. If the same form of transmission is used 
as for the S;vept Frequency Mode (i.e. three sets of three pulses, prf of 300 pps), 
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lu frequencies would Ik* cycle through in 0.3 seconds giving a resolution along the 
orbit of a few kilometers. This mode will In* particularly useful for studying var- 
ious types of ionospheric structure such as TIO's, patches of spread F distant 
from Spaeelab and bubble structures of small transverse dimensions. 

Figure *< shows a Doppler coded group delay vb. time record produced by an 
irregularity displaced from the orbital path. At closest approach, the Doppler 
shitt is zero and the group dcluy a minimum. 

I)c|M*nding on the number of frequencies being cycled inrough, and their spac- 
ing. it may be possible to invert the vertical incidence data and determine the 
electron density contours of the background ionosphere. Kay tracing could then 
Ik* used to determine the location of the irregularity shown in Figure 6. 

Vlternatively, the location may be estimated as follows. If the earth's mag- 
netic field is neglected, the group range of an oblique echo may be related to the 
group height at vertical incidence, by (Calvert and Cohen, lJHil), 

hj, 2 « h' v 2 ♦ R 2 (8.1) 

l\/ f v B h t /f s < 8 - 2 > 

where It is the horizontal range to the irregularity, h ^ is the apparent range at 

the frequency l' H , and h' is the apparent vertical range at f , the frequency re- 
flected vertically at the same height as the oblique ray at frequency f R , (see 
Figure 7). Now according to Mart.vn's Theorem 

f v = f R cos 6 (8.3) 

where is the angle of Incidence at the m ■ 1 level. 

For a vehicle imbedded in the ionosphere, Snell's Law gives 

co » f ... B * in<? o < 8 * 4 > 

If the vehicle is moving horizontally with velocity v, at an angle of y to the 
plane containing the ray path (Fig. 8), the Doppler shift may be written as 
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. 2f . « 

{ e T^..« v ,,n6 *..t 


cot 


(H.5) 


Using (0.1), (8.2) and (8.3) Hives 

h| iind ( ■ R (8.6) 

Using (8.6), (8.4) and (7.5) we obtain 


K cos f 


r 2iil 1 h 


(8.7) 


Hence, measurements ot' apparent range ligand Doppler shift *. f enable the pro- 
jection of the irregularity true range along the vehicle orbit to be obtained. Con- 
sequently, the location of irregularities can be estimated from portions of Fixed 
Frequency records, or even Swept Frequency records. Thai is, it is not neces- 
sary to observe either a crossing of zero Doppler shift or a minimum In apparent 
range in order to determine the position of the irregularity. 

Doppler Spectra. This mode requires operation at one frequency for ~0.6 
seconds. It is consequently a relatively time-consuming mode of operation. It 
is best used by alternating with swept frequency modes or interspersing it at 
fixed time intervals in the FixeJ Frequency Mode. Five or ten selected frequencies 
could be used taking approximately a total time of 2.5 or 5 seconds to obtain the 
Doppler Spectra. 

A future development to consider is one in which two receivers are used, 
one for Doppler spectra and the other for the Swept Frequency and Fixed Fre- 
quency Modes. A second transmitter would not be needed, since the transmitter 
would only need to be capable of transmitting two interlaced sequences of pulses 
with a sufficient frequency difference so that each receiver would only be sensi- 
tive to one or other of the sequences. 
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Control Program. The operation mode of the instrument and selection ol 


operating parameters (frequency, p.r.f., step size, etc.) will be under computer 
control. The operation modes will be selected or changed by providing the ap- 
propriate input parameters to the control program. The appropriate input for 
several types of operation sequences will be available from the computer memory. 
The payload specialist will be able to enter other operation sequences. It is rea- 
sonubly straightforward to program several modes of operation and allow a wide 
choice of parameters (e.g. sequence of modes, swept frequency range, frequency 
step, number of fixed frequencies, etc.) and this approach has been used success- 
fully with ground based ionosondes (e.g. Devlin, 11)78). 

if. Data Rates 

In the discussion, thus far, all techniques have been based on measurements 
ol the in-phase and quadrature amplitude compoents of the signal. In section 3, 
a lo kHz bandwidth was assumed so that the in-phase and quadrature amplitudes 
need lo be sampled at a 20 kHz rate. Each amplitude is assigned either an 8-bit 
or 12 bit word, so that the measurements require a data rate of 40 k words/sec. 

The maximum p.r.f. suggested here is 300 p.p.s. and for each pulse the 
frequency and AGC level must be recorded. The AGC requires one word or less 
(see Section 2). In the MF anil II F range, frequency steps as small as 20 kHz 
have been used in this discussion. It is unlikely that other modes of operation 
would require a frequency stop of less than 1 kHz. If the frequency is to be 
specified to 1 kllz, we can divide each range of 1 megahertz into 4 sub-ranges 
of 230 kHz. The basic range of 25o kHz can be specified by 8 bits, and to de- 
termine the correct sub-range over 30 MHz requires 7 bits. Thus two words 
are required to specify the frequency. 
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Hence, the overall word rate 1 b 40.9 k words/sec. A rate of 41 k words/sec. 
would allow loo worda/aec for other information, such as universal time, and 
operating mode to be specified periodically. For s bit words, this corresponds 
to a bit rate of 3.28 x 10 s bits/ sec. 

10. Oblique Measurements 

The types of measurements discussed previously can all be done by trans- 
mitting from one vehicle and receiving on another. Some advantages of doing 
this have been discussed by Calvert (1976) and will not be repeated here. Only 
one aspect of combining Doppler ;ind time delay measurements between two 
vehicles will be discussed. 

From measurements on a single vehicle, the apparent range and angle of 
arrival (as deduced from the Doppler shift) cun be determined. The apparent 
range, R' defines a sphere in R' space on which the echo lies. Since echoes re- 
turning from in front of the vehicle have a Doppler shift of different sign to those 
returning from behind the vehicle, the angle of arrival defines a circle on the 
sphere on which the echoing region must lie (Figure 9). In the absence of further 
measurements, more precise location of the echoes requires assumptions to be 
made and in some instances these may be reasonable. For example, if the angle 
of arrival corresponds to propagation along the magnetic field direction, it would 
be entirely reasonable to assume that propagation was along the field line rather 
than at an angle to it. In other situations, however, it is likely that simplifying 
assumptions could not be made with confidence. 

If measurements are made between two vehicles, then the apparent time 
delay, which is equivalent to an apparent path length, P' defines an ellipsoid 
upon which the echoing region must lie (Figure 10). Since the echoing region can 
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be confined to a circle by the single vehicle measurements, then, in the case of 
|K»int-like echoing regions, the echoing region must lie on the intersection of the 
circle with the ellipsoid. The cross section of the ellipsoid in the plane contain- 
ing the circle will be an ellipse. The ellipse and the circle can have from 2 to 4 
intersections (Figure 11). When there are 3 or 4 intersections, some of the 
ambiguity may be reduced by using the Doppler shift information for propagation 
between the two satellites as follows. Since the apparent range from the vehicle 
containing the monostatic measurements is known, the apparent range from the 
other vehicle is jus? the difference between the monostatic range and the apparent 
path length between the vehicles. In this way the location of the echoing region is 
either determined uniquely or known to be one of two possible locations. With 
the geometry of the situation now known, the expected Doppler shift between the 
two vehicles can be calculated and the degree of agreement with the measured 
Doppler shift will indicate the extent to which the reflecting region acts like a 
point reflector. If there is significant disagreement, some knowledge of the ex- 
tent of the reflecting surface may be estimated by considering different models. 

11. Measurements of the Horizontal Wave Vector of Gravity Waves 
Combining Monostatic II. F and Bistatic VHF Measurements 
Suppose Spaeelab contains the H.F. monostatic system and that the sub- 
satellite has a VHF receiver to measure the Doppler of a VHF signal trans- 
mitted from Spaeelab. For simplicity it is assumed that the two vehicles are 
traveling northward and separated by a distance L in the east-west direction 
(Fig. 12). Both vehicles travel horizontally and with the same velocity. 
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Suppose that there is a wave-like perturbation in the electron density and 
that in the horizontal plane containing the space vehicles the variation in elec- 
tron density can be written as 

-f) <“- 1 > 

l.e. wavelike perturbation with apparent wavelengths in the x direction and 
in the y direction (Fig. 12). Fronts of constant phase occur where y "(k y A»)x 
i.e. phase fronts are tilted to the x axis at an angle 6 * artan K y /k K and the wave 
vector is at an angle v/2 - 6 to the x axis and the true wavelength \ T is given by 
lAj ■ l/\ 2 • 1/A 2 . Spacelab moves along the line x • 0. 

The variation in N e along the orbit of Spacelab is of the form 


AN( x . y ) = AN 0 c 


os 2rr 


AN 




= AN 0 cos2n 


= A N 0 cos 2-r 



( 11 . 2 ) 


and so \ y can be determined. 

The sub-satellite moves along the line x = L. Consider the phase path P, 
of a VHF ray propagating from Spacelab to the sub-satellite. 


P = 





cos a ds 


For VHF we may consider ray paths to be line of sight, so that 

P = f M dx 

Jo 

Ignoring the magnetic field (again an excellent approximation at VHF), 

M 2 = 1 -X 
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u ~ 1 - — Xsincc*X<< 1 
2 


= 1 - — K — - sinceX - — - 
2 f j f j 


Thus 

J * L 

N dx 
o 

dP_ 1 K d f L ., 

dt ‘ 2 f 2 dt J o N * 

where K is a constant. The Doppler shift is given by 


c dt 


,-i JLA f n 

2 f>c d *J 0 * 


(Note that dL/dt has been assumed to be zero, i.e. we have assumed that the two 
vehicles maintain their spacing with time. II this is not true then 


dP_ dL _l JK _d C\ 
dt dt 2 f 2 dt J o * 


However if two frequencies are used, as in the proposed ixperiment (e.g. Grossi 


and Gay, 1975) we have 


jl a r N d, 

dt dt 2 r 2 dt J e 
1 1 •'0 

i.i‘I ± f L n * 

dt dt 2 fi dt J ' 
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and 


^ dP 2 1 „ 
dt " dt "2 


h hd U [ l n 

v / d ' J. 


dx 


Thus the differential Doppler shift can be used to eliminate dL/dt. Even with 
measurements on only one frequency, dL/dt is likely to be a slow variation com- 
pared with the effect of a gravity wave and so could be filtered out). Thus we 
are concerned with determining 




dx 


Now N e ■ N h . C |, fround * AN(x,y). Consider a particular instant when Spacelab is 
at (o,y) and sub-satellite is at (x.y). The instantaneous value of J N,dx will be 
given by 


N ' = f ( N 


b + AN 0 cos 2 


•fe-i) 


dx 


Now y ■ vt, so N t will vary with time according to 

n t = N - • AN « { ,in 2 ”(r7 - 77 ) - * in (-X 7 )} 


and 


‘* N t am / « / L vt\ 0 / vt \ 

— — = -V — AN ft < cos 2 tt[ -cos 2 n — 

dt \ y 0 \ U* S / \ *y / 


(11.3) 


which may be written as 


dN 

"d 


— = 2v~AN 0 sin2t7 s i n (11.4) 

t K y 0 \ 2 KyJ K x 
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Comparing (11.2) and (11.4) we see that the VHF Doppler variation will have the 
same period as AN n but there is a pha..e shift of -LA,. Thus if l, is known, 
we can determine K t and hence \ T . (Note however there will be an ambiguity 
in k, i.e. two possible directions differing by 180°). 

From (11.4) we see that if \ y -• «, there is no ilN T /dt due to the wave. This 
is the expected result. When \* - ® , the phase shift approaches 90° as we ex- 
pect and since 


sin(77L/\x) (TrL/k*)-*! 


£t^vL 

(It \y 


AN 0 s in 2n 



i.e. dN,/dt depends only on \ y as we expect. Suppose a monostatic HF system 
or a direct measurement probe is on the sub-satellite. It would measure vari- 
ations given by 

AN = AN 0 cos 2 v ( — - — 

i.e. a phase shift of 2v L/\ x compared with Spacelab measurement. However, 
there are ambiguities in X* such that 


2 ttL /A. x = 2iit 7 , n = 1 , 2 , - - - 
i.e., L/X x = n \ x 

In the case of the VHF Doppler experiment there are ambiguities in phase when 

v L/X x = 2nv 
i . e., L = 2nX x 

However, the amplitude term is proportional to \ x and so the ambiguities c m 
be removed. 
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Note: The Doppler experiment hue zero amplitude when L ■ n\ bo that the best 
results would be obtained by combining the VHF Doppler with a local measure- 
ment (either by H.F. system or probe, «n« sub-satellite. 
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Figure 1. Reflected ray path in a plane 
stratified ionosphere 



Figure 2. Reflected ray paths from an 
irregularity embeddud in a plane strati- 
fied ionosphere 
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Figure 3. Phase shift due to Doppler effect, A/, between 
two echoes of two pulses transmitted at the same fre- 
quency but at a time interval of 1/300 sec. 
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Figure 4(a). Possible echo paths when 
Spacclab is within an equatorial bubble 
or duct and when a field-aligned irreg- 
ularity patch exists below Spacelab 
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Figure 4(b). Ionogram obtained for the situation 
shown in Figure 4(a) 
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Figure 5. Doppler coded ionogrutr. with a range window of 500 km 
obtained for the situation shown in Figure 4 
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Figure 6. Doppler coded P' vs. t signature 
produced by an irregularity embedded in 
the ionosphere below Spacelab 
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Figure 7. Flay paths having related delays and frequencies 
as given by Calvert and Cohen (1961) 



Figure 8. Ray path geometry for an irregularity 
displaced from the orbital plane 
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Figure U. Echo surface determined from monostatlc group 
delay and angle if arrival (Doppler) measurements 


p‘ = p; + p ‘ 



Figure 10. Echo surface defined by the group delay 
from Spacelab to a sub-satellite 
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Figure 11. Intersections of monostutic and 
bistatic echo surfaces 



Figure 12. Geometry for Spacelab and a sub-satellite 
moving northward through a region containing a grav- 
ity wave 
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